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We study numerically an atomistic model which is shown to exhibit a one-step crystal-to- 
amorphous transition upon decompression. The amorphous phase cannot be distinguished from 
the one obtained by quenching from the melt. For a perfectly crystalline starting sample, the transi- 
tion occurs at a pressure at which a shear phonon mode destabilizes, and triggers a cascade process 
leading to the amorphous state. When defects are present, the nucleation barrier is greatly reduced 
and the transformation occurs very close to the extrapolation of the melting line to low tempera- 
tures. In this last case, the transition is not anticipated by the softening of any phonon mode. Our 
observations reconcile different claims in the literature about the underlying mechanism of pressure 
amorphization. 

PACS numbers: 



I. INTRODUCTION 

Since the first reports on pressure-induced amorphiza- 
tion (PIA), there has been an increasing interest in the 
phenomenon that has been observed in a large class of 
systems, both upon compression and decompression (see 
Ref. [1] for a review and references on particular ma- 
terials). The experimental interest stems from the fact 
that PIA occurs in some extremely widespread and im- 
portant materials, namely water— and quartz and its 
polymorphs-, and also because PIA provides a novel 
route to the synthesis of amorphous materials, in addition 
to the traditional technique of quenching from the melt. 
^From a basic point of view, there are controversies about 
the mechanism of PIA and the nature of the amorphous 
phase, in particular regarding whether it is equivalent or 
not to the one obtained by quenching^*. Since in most 
cases crystalline phases that undergo PIA show reen- 
trant melting in the P-T plane, it has been suggested* 
that PIA is nothing but melting at temperatures below 
the glass temperature of the supercooled fluid. Other 
studies (specially numerical ones 6 ' 7 ' 8 ' 9 ) have instead em- 
phasized the relationship between PIA and mechani- 
cal instabilities. It has been observed in fact, that in 
many cases PIA is triggered by the softening of a shear 
phonon modeiS. In addition, some memory effects^* and 
anisotropic properties*, although controversial^*, show 
that many times what appears to be an amorphous phase 
preserves within its structure signatures of the parent 
crystalline phase. 

In view of the broad phenomenology briefly stated 
above, it is highly desirable to take advantage of model 
systems in which amorphization can be studied in a 
transparent way, allowing to look in detail into the mech- 
anisms and characteristics of the transformation. This 
motivates the present work, in which we study a simple 
two-dimensional system of identical point -like particles 
interacting through a specially devised two-body poten- 



tial. The simplicity of the model allows us to study large 
systems with some amount of defects, and observe di- 
rectly the crucial role they play in the transformation. 

In the model we study, there are a few different crys- 
talline ground states depending on the applied pressure. 
We present here the results of the evolution of the most 
compact structure (stable at the highest pressures) upon 
pressure releasei*. As we will see, PIA at temperatures 
in which particle diffusion is negligible is always related 
to mechanical instabilities. For perfect lattices, PIA is 
reflected in the softening of a shear phonon mode. This 
leads to local distortions that produce the destabiliza- 
tion of new vibrational modes, leading to a cascade of 
instabilities 8 , which drives the system toward an amor- 
phous structure. However, in the presence of defects, 
localized vibrational modes exist, that may become un- 
stable before any extended vibrational mode does. This 
favors the nucleation of the amorphous phase at pressures 
much closer to the thermodynamic equilibrium value be- 
tween the crystalline and amorphous phases. 

The work is divided as follows. In Section II we present 
the model. The results at zero temperature are contained 
in Section III. In Section IV we give evidence that the 
disordered samples obtained can be called trully amor- 
phous. In Section V we present results at finite tem- 
peratures, and Section VI contains some discussion and 
conclusions. 



II. THE MODEL 

An isotropic, purely repulsive interparticle potential is 
used, with a strict hard core at a distance ro plus an 
almost triangular repulsive shoulder. The pairwise inter- 
action potential V(r) between two particles separated a 
distance r is given byi* 

V(r) = oo for r < ro 
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where r\ — 1.17315ro, and ro sets the length unit. The 
potential is plotted in the inset of Fig. ^ This kind 
of potential has been previously used to systematize 
the anomalous properties of tetrahedrally coordinated 
materialsi2ii£*i&. Then we expect they are also appro- 
priate to study amorphization under pressure, since this 
phenomenon occurs for most of these materials. Its cru- 
cial characteristic is the existence of two possible equilib- 
rium distances between particles. 

The system is simulated by standard molecular dynam- 
ics in the NVT ensemble with periodic boundary condi- 
tions. At a given volume, the quantities of interest are 
evaluated, and the volume is changed in steps of the or- 
der of 0.01% by rescaling all coordinates of the particles 
and size of the simulation box. Temperature is fixed by 
rescaling the velocities of the particles whenever neces- 
sary. As the amorphization process implies the existence 
of mechanical instabilities, we typically observe that ki- 
netic energy tends to increase during amorphization. In 
a real situation this energy transforms into heat. Here 
we simply eliminate it by the mentioned rescaling pro- 
cedure. Pressure is calculated by a direct evaluation in 
terms of the interparticle forces. We choose to model 
the system at constant volume in order to survey all re- 
gions of the volume-pressure curve, including those that 
would be unstable in constant pressure simulations. The 
results to be presented correspond to a two-dimensional 
system, to facilitate a direct visualization of the parti- 
cle configurations. We should mention however, that the 
same phenomenology was observed in three-dimensional 
samples with the same interparticle potential. 



III. ZERO TEMPERATURE RESULTS 

We take as the starting configuration the one corre- 
sponding to the triangular lattice with lattice parameter 
~ r , which is stable at high pressures. In Fig. Qwe 
present the results of simulations at T = in a system 
of 2800 particles. The thin continuous line corresponds 
to the v-P relation assuming the system remains always 
triangular. The reentrance of this line is a consequence 
of the particular form of the interaction potential. The 
numerical results for a perfect lattice (full circles) follow 
this line up to some maximum volume, at which they 
abruptly depart from it. This is the pressure in which a 
phonon with vanishing energy first appears in the trian- 
gular structure (see Fig. [2J . Analytical evaluation shows 
that this instability occurs when 3V'(a)/a + V"(a) = 0, 
where primes denotes derivatives of the potential, and a 
is the lattice parameter. This expression is valid as long 
as the interaction between next nearest neighbors is zero, 
as it is in the present case. The critical volume v c and 




FIG. 1: Evolution of pressure P upon increasing of the spe- 
cific volume v, at zero temperature. The thin continuous line 
is the expected evolution if the system remains always trian- 
gular. The dotted lines marks the ideal stability limit (v c , P c ), 
at which a phonon energy vanishes. Dots indicate results of 
simulation for a perfect lattice (full circles), a lattice with 
a single vacancy (stars) and a system with grain boundaries 
(open circles, see Fig. Letters correspond to snapshots in 
Fig. [|] The inset shows the interparticle potential V(r). 



pressure P§ are indicated by dotted lines in Fig. ^ an d 
they are fully compatible with the numerical results. The 
unstable phonons turn out to be shear phonons with the 
wave vector k oriented perpendicularly to one of the three 
most compact directions in the lattice. Note that for the 
interparticle potential we use, the longitudinal phonon 
branch along these directions is of the type ~ sin(fc), and 
all shear phonons -irrespective of the absolute value of 
k— become zero energy at the same point. Then, in the 
present case the instability pressure is macroscopically 
signaled by the vanishing of the stress dependent shear 
modulus of the materialize, fi = (C u - C 2 2)/2 - P = 0. 

At v c there is a sharp and abrupt increase of pressure^. 
Snapshots of the particle configurations give clues of what 
happens in the system. In Fig. [«5fa) we see the triangu- 
lar configuration just before the instability. Figure E[b) 
shows that after changing the volume a very small quan- 
tity, an instability has propagated in the system. We 
stress that the evolution of the system from the config- 
uration in Fig. GJa) to that in Etb) is triggered by a 
very small volume change, and is just the steepest de- 
scend evolution toward a local energy minimum of the 
energy landscape. In this evolution, pressure recovers to 
a higher value indicated as P* (dashed line) in Fig. ^ 
This value is roughly maintained2£ upon further volume 
increase up to v/tq ~ 1.27, which is close to the volume 
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FIG. 2: The dispersion relation of phonons along one of the 
directions perpendicular to the densest planes of the trian- 
gular structure close to the instability pressure. Continuous 
lines correspond to the repulsive potential Eq. Q , and dashed 
lines to the same potential plus an attractive term (Eq. @) 
that includes second nearest neighbors interactions. Curves 
correspond, from top to bottom, to \P — P§\/{r^/eo) = 0.01, 
0.001 and 0.0001. Note that since the repulsive interaction 
reaches only the first neighbors, the full branch has an ana- 
lytical form of the type ~ sin(fc), and all phonon vanish at 
the same pressure P = P§. For the case with second nearest 
neighbors interaction the instability occurs at the zone border 
phonon ka = 27r/\/3. 




at which the disordered regions have taken over the whole 
system. If volume is increased further, pressure decreases 
noticeably. 

The thermodynamic equilibrium pressure P e at T = 
between the crystalline and disordered structures is in- 
dicated by the arrow in Fig. ^ The value of P e was 
calculated by enthalpy evaluations of the triangular lat- 
tice and the amorphous structure obtained by quenching 
from the melt. Thus, P e represents the natural extrap- 
olation of the melting line to zero temperature. We see 
that P e is close to (though a bit higher than) P*. 

It must be noticed that constant pressure simulations 
would have shown at Pq an abrupt transition between a 
completely ordered and a completely disordered system, 
while P* would have remained hidden. However, P* has 
a clear physical meaning as the transition pressure once 
the disordered phase has been nucleated. This suggests 
that if nucleation centers are present in the system, P* 
will be experimentally accessible as the actual transition 
pressure. In fact, Fig. ^ shows also the v—P evolution 
for a system with a single vacancy and for a system with 
grain boundaries. In the former case, the instability oc- 
curs before v c is reached, and in the latter case it even 
occurs without any pressure reentrance^i. Note that the 
three curves tend to coincide after the first destabiliza- 
tion of the original lattice. Two snapshots of the system 
for the poly-crystalline case [Fig. |3Je) and (f)] show how 
the disordered phase grows from grain boundaries, that 
provide nucleation centers for the transformation. These 



FIG. 3: Snapshots of the systems at the points indicated in 
Fig. Q (only about a quarter of the full simulated system is 
shown). Upper four panels correspond to a perfect crystal 
starting sample, whereas the two lower ones are from a poly- 
crystalline sample with grain boundaries and vacancies. 



results are a clear indication that defects are very effec- 
tive in lowering the energy barrier for the transition, and 
also explain why, in experiments, transformation pres- 
sures close to the thermodynamic values are usually ob- 
served. The lowering of the energy barrier due to defects 
and its relation to thermodynamic melting has been dis- 
cussed by Mizushima et alm& in the context of crystal- 
crystal pressure-induced transformations. Note that in 
the case in which the transformation is triggered by de- 
fects, the phonon spectrum of the system (and in par- 
ticular the elastic constants) gives no indication of the 
instability that is about to occur. This is particularly 
obvious in the case of a single vacancy: a single defect 
cannot modify the phonon spectrum of an infinite sample, 
and then its effect it is not seen in the elastic constants, 
but it produces a finite change in the value of the critical 
pressure Pq— . 

It was already mentioned that for the interaction po- 
tential used, the whole branch of shear phonons become 
zero energy at Pq. We want to emphasize however that 
amorphization it is not related to this degeneracy, as it 
can also be triggered by a single phonon becoming un- 
stable. In fact, we did simulations with a modified po- 
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FIG. 4: The structure factor of the amorphous obtained by 
decompression of a perfect crystalline sample and that ob- 
tained by a quench from the liquid phase, at the same volume 
v/ro = 1.289. They coincide within the numerical precision. 



tential in which an attractive term V a reaching up to 
second nearest neighbors was included. The actual po- 
tential used is the one of the previous section plus an 
attractive term given by 



Va(r) 
Va(r) 
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The instability corresponds now to the vanishing of the 
energy of edge zone shear phonons only, as indicated in 
Fig. [3 The amorphization of the system occurs precisely 
at the point where this phonon becomes zero energy, and 
it was observed to be similar to the previous case where 
the full shear phonon branch becomes unstable at the 
same pressure. This rules out an amorphization mecha- 
nism in which at the instability, arbitrary combinations of 
the unstable phonons generate disordered movements of 
the atoms, and favors a cascade mechanism as described 
in the discussion section below. 



IV. COMPARISON BETWEEN 
DECOMPRESSED AND QUENCHED 
AMORPHOUS 

The claim that the disordered structure we obtain 
upon decompression is really amorphous is supported by 
the following facts. We find no evidence of orientational 
order in the structure factor <S(k), even in the case in 
which we start with a mono-crystalline sample. More- 
over, a quantitative comparison with the structure factor 
of a system quenched from the liquid (prototypic of an 
amorphous sample) shows that they are indistinguishable 
(Fig. |3} . In addition, we have failed to find any single 
systematic difference between the two amorphous. 



V. FINITE TEMPERATURE EFFECTS 

If temperature is not strictly zero, even the perfect 
lattice amorphizes before the ultimate mechanical insta- 
bility limit is reached. Simulations at finite temperature 
with perfect samples show that the amorphization occurs 
at a temperature dependent critical pressure P°(T) (such 
that P C (T = 0) = Pq). The temperature dependence of 
the critical pressure originates in the fact that a ther- 
mally activated process may destabilize the soft phonon 
before it actually becomes zero energy. 

It is instructive to see in some more detail how 
this destabilization occurs. Consider the normal modes 
(phonons) of the system. We will study the case in which 
only a single phonon with wave vector k perpendicular to 
one of the high density planes has non-zero amplitude. 
We expand the energy of the lattice in powers of the 
amplitude of the phonon oscillation, going one order be- 
yond the harmonic approximation. Let us call Si and St 
the longitudinal and transversal oscillation amplitude of 
that phonon. As we will see there is a nontrivial coupling 
between the two at the instability. 

Taking into account that our potential produces in- 
teractions only among first neighbors in the lattice, the 
energy per particle eo of the system up to third order 
in the amplitudes can be written after some lengthy but 
straightforward calculation as 



eo 




(3) 
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The first line in J3J| is the harmonic contribution, the 
second is a cubic term in the displacements (there is an 
additional cubic term proportional to Sf , but this can be 
shown to be not relevant for the analysis below). The 
system becomes unstable when the coefficient of S\ van- 
ishes, providing again the condition 3V'/a + V" = 0, 
which defines the T = values Pq and v c . Close enough 
to this instability point, 3V'/a + V" is proportional to 
P — Pq an d to lowest order all other coefficients can be 
considered to be constants. Then generically, the energy 
can be written as 
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The existence of a Si Si term is very important, as it indi- 
cates that for P slightly larger than Pq there is a saddle 
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point (actually two, one with positive and the other with 
negative S t ) close to S t = 0, Si = 0, than can be deter- 
mined requiring stationarity of The result is 
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where the upperscript indicates the values at the saddle. 
If the system reaches this saddle it can escape from the 
local minimum at the origin, namely, the system desta- 
bilizes. The energy barrier h for this process is obtained 
by reinserting into Q. The result is: 



h = 



4A 2 B(P - P c 
D 2 



(0) 



Note that the barrier is exactly the same whatever the 
wave vector of the phonon considered (this is no longer 
true if further neighbors interactions are included, but 
the possibility to escape through the jump of a barrier 
remains). The instability mechanism is then driven by 
the Sf6i term in the energy, what implies a not trivial 
coupling of transverse modes (the one actually having 
vanishing frequency) and longitudinal ones. At finite 
temperatures the barrier h can be surmounted. Nuclc- 
ation theory tells that escape time r is proportional to 
exp (h/T). Assuming the preexponential factor is a con- 
stant To (this is certainly not true, but dependences on r 
become weak in the final result) we obtain that if a time 
to is given, the system will overcame the energy barrier 
if to r, and from here we obtain the formula for the 
temperature necessary to escape a given barrier, namely: 



T> hln- 1 (t /t ). 



(7) 



For practical purposes the logarithmic factor can be usu- 
ally taken to be approximately 0.1. As h ~ (P — P§) 2 , 
we obtain that the critical pressure P C (T) increases as 
T 1 / 2 at finite temperatures, a behavior that is well repro- 
duced in the simulations (Fig. [5J). The P C (T) line is then 
the pressure amorphization line for a sample without de- 
fects, and can be properly interpreted as the 'mechanical 
melting' line of this system. In the absence of defects it 
marks the limit on which the crystalline phase destabi- 
lizes (purely mechanically at T = 0, and by a thermally 
activated process at T ^ 0). Note however, that the 
elastic constants of the system, in particular the shear 
modulus fi, do not extrapolate to zero at P C (T), except 
at T = (see Fig. |HJ|. In fact, the numerical evaluation 
by means of fluctuation formulasiSiSi shows that the ef- 
fect of temperature in the elastic constants is very small. 
At T = the stress dependent shear modulus \x vanishes 
at Pq , as we already knew from the phonon dispersion 
relation, but at T ^ extrapolation is consistent with 
a vanishing of \i very close to Pq , and certainly not at 
P C (T). The thermally activated process that at finite 
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FIG. 5: Evolution of P c (full circles) and P* (full squares) 
as a function of temperature (bars in P* indicate the whole 
range in which transformation occurs). The inset shows in 
a log-log plot that AP C = P C (T) - P C (T = 0) increases as 
T 1//2 at finite temperatures. The extrapolation of the melting 
line (from calculations at higher temperatures, two points are 
seen) is indicated by the dashed line. It goes to P e at T — 0. 



temperatures leads to the instability at the higher pres- 
sure P C (T) is a rare event that it is enough to occur once 
to completely destabilize the system. This process is not 
captured in the value of the elastic constants. 

Up to here the results for samples without defects. For 
defective samples the amorphization line P*(T) is only 
weakly dependent on T. In addition, it is very close to 
the extrapolation of the melting line obtained from sim- 
ulations at higher temperatures (see Fig. EJ. P*(T) and 
P C (T) are well different from each other, particularly at 
very low temperatures. Defects are responsible for this 
difference, and provide the link between the 'mechanical' 
and 'thermodynamical' scenarios for pressure amorphiza- 
tion. 



VI. DISCUSSION AND CONCLUSIONS 

The experimental observation that in some materi- 
als PIA occurs roughly along the extrapolated melting 
line has motivated the suggestion^ that 'thermodynamic 
melting' is the underlying driving force for amorphiza- 
tion. Our results show that amorphization at zero tem- 
perature is always a phenomenon related to mechanical 
instabilities. In perfect lattices, mechanical instabilities 
reflect necessarily in the phonon spectrum of the mate- 
rial, and amorphization can be easily related to a me- 
chanical process. However, for defective lattices, the me- 
chanical instability that triggers amorphization is asso- 
ciated with localized vibrational modes around defects, 
and this instability occurs typically without any notice- 
ably signature in the phonon spectrum, i.e., the mechan- 
ical nature of the process is more subtle. We emphasize 
that a single defect is sufficient to trigger the amorphiza- 
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FIG. 6: Variation of the stress dependent shear modulus fi 
with pressure at different temperatures. For T = 0, fj, van- 
ishes at the critical pressure P C (T = 0), but for increasing 
temperatures fi does not tend to zero at the corresponding 
critical pressure P C (T), showing that the transition it is not 
macroscopically captured by the elastic constants. The nu- 
merical values of P C (T) are indicated by the vertical arrows. 



tion of a whole macroscopic sample^. 

How this first instability is able to transform the highly 
ordered original lattice into a fully amorphous configura- 
tion is not easy to understand, specially considering that 
this evolution is fully deterministic from a mechanical 
point of view, and can be described as the steepest de- 
scendant path of the configuration point of the system 
onto its energy landscape. A deep understanding of this 
process will surely shed light also onto the very defini- 
tion of what an amorphous material is, typically char- 
acterized by what it lacks, instead of what it possesses. 
Here we make only the following considerations?. The 
first instability of the lattice (at T = 0) can always be 
described as the vanishing of the frequency of one of the 
normal modes of the system. This normal mode may be 
a phonon for a perfect lattice, or a localized mode for 
a defective lattice. When the coordinate characterizing 
the destabilizing mode grows, it typically goes beyond 
the applicability limit of the linear theory, and the whole 
mechanical stability of the system has to be reanalyzed. 
Then, although the analysis of Section is appropriate 
to determine the barrier for the system to scape from its 
metastable equilibrium, it cannot be continued when the 
system actually overcomes the saddle without including 
the interaction with other modes. For instance, it occurs 
in our system that once the coordinate of the unstable 
phonon starts to increase, it couples to other phonons 
and the results is that a shuffle of two pieces of the ma- 
terial develops. In this way, the instability generated by 
an extended object (the unstable phonon) produced a lo- 
calized (one-dimensional) defect on the latticed. But it 
turns out that, in the present case, the new configuration 
of the lattice is not stable either, and a new distortion 
occurs spatially close to the region already distorted (see 
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FIG. 7: The first steps in the cascade process from a per- 
fect to an amorphous lattice: The original perfect lattice be- 
comes unstable when the energy of a phonon vanishes, and 
this drives a transition to a sample with a shuffle (a). This 
configuration however is not stable and spontaneously evolves 
toward the configuration in (b). But this is not stable yet, and 
the process continues until equilibrium is reached when a fi- 
nite fraction of the system is transformed (compare with Fig. 
(2) (a) and (b)). 



Fig. UJ . In the present case of amorphization, it is likely 
that a cascade of these processes occurs that leads to an 
amorphous structure. The spatially localized nature of 
the process reflects in the fact that for partially trans- 
formed samples, transformed and untransformed regions 
are spatially separated, as seen in Fig. [31 We however 
emphasize that another possibility could be that after the 
first phonon instability and the first shuffle, the system 
reaches a mechanically stable configuration, then further 
increases of the volume will produce independent shuf- 
fles in the samples, along the three directions related by 
symmetry. If this is the case, then this single phonon 
will drive a transformation (which has the character of a 
martensitic transformation) to a new (poly-)crystallinc 
phase. This and other intermediate cases, that in our 
model occur for different parameters, will be discussed 
separately. 

In summary, we have arrived at a picture of the mecha- 
nism of PIA, reconciling mechanical and thermodynam- 
ical melting in a unifying scheme. The present model 
exhibits a one-step crystal-to-amorphous transition to a 
truly amorphous state that does not differ from the dis- 
ordered structure obtained by quenching from the melt. 
Once the transformation is triggered, it is followed by a 
cascade process that produces the amorphous structure. 
For a perfect lattice, and vanishingly small temperatures, 
the transformation takes place through a shear phonon 
instability at a pressure far apart the thermodynamical 
equilibrium pressure. For defective lattices, the mechani- 
cal instability shows up by the destabilization of localized 
vibrational modes, close to the thermodynamical equilib- 
rium pressure. In cither case, the transformation is al- 
ways triggered by mechanical instabilities. Temperature 
can produce the transition to occur before the mechanical 
stability limit is reached, due to a process of thermally ac- 
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tivated jump over a barrier. In real samples with appro- 
priate defects, experimental transition pressures will be 
seen close to the extrapolation of the melting line. These 
facts reconcile different claims in the literature regarding 
the nature of the PIA of materials. We remark that, al- 
though we presented results for a two-dimensional model 
system, preliminary results in three-dimensional samples 
qualitatively agree with the phenomenology stated above, 
and will be presented elsewhere. We finally stress that 
these considerations have potential applications in the 



study of melting mechanisms of crystalline materials 26 . 
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